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Volumetric flow-rate (L/hr) 

Haldane kinetic parameters 
No. of stages 

Ratio of two successive volumes of the reactor 
Substrate concentration (g/L) 

Feed substrate concentration (g/L) 

Initial substrate concentration (g/L) 

The substrate concentration at which the reaction rate 
is a maximum 
Time (hr) 

Final time (hr) 

Duration of sub- interval (hr) 

Biomass concentration immediately after the addition of 
each pulse 

Volume of the reactor (L) 

Initial volume of the reactor (L) 

Final volume of the reactor (L) 

Volume of pulse (L) 

Biomass concentration (g/L) 

Initial biomass concentration (g/L) 

Final biomass concentration (g/L) 

Constant yield 


Greek: 

li : Specific growth rate (hr 1 ) 

I n : Specific growth rate constant 



Subscripts: 


0 : • Initial 

f : Final 

F : Feed 

. th' . 

1 : 1 stage 

Superscripts: 

: At the end of each sub- interval 

+ : After the addition of the pulse, or at the begining of a 

sub interval . 

* : maximum 



Abstract 


In this work we analyse the optimal control problem of a 
biochemical reactor sustaining a reaction characterised by 
substrate inhibition. The control variable is the substrate feed 
rate. We approximate the feed rate policy by discrete pulses. 
The entire period of operation is divided into sub- intervals . The 
substrate is assumed to be added at the of each sub- interval 
as pulses. 

The optimisation problem is solved analytically when the 
initial condition is stoichiometric. A sequential quadratic 
programming is used to determine the optimal policy when the 
initial condition is non- stoichiometric . We consider two cases i) 
when all sub-intervals are of equal duration ii) when the optimum 
duration of each sub-interval has to be found. The effect of 
initial conditions on the optimal profile is also investigated. 
We consider four possible modes of the operation of the reactor 
when (i) the final volume is fixed, final time is fixed, (ii) the 
final volume is free, final time is fixed, (iii) the final volume 
is fixed, final time is free, (iv) the final volume is free, final 
time is free. 





CHAPTER 1 


INTRODUCTION 

Bio- chemical reactions characterised by substrate inhibition 
exhibit an optimum performance when the reactor is operated in the 
fed-batch mode. Here the substrate is fed to the reactor over a 
period of time to ensure that the reaction rate is close to the 
maximum in the reactor (Cazzador (1988) , San and Stephanopolous 
(1984) ) . 

Pontryagin's (1977)' maximum principle has been used to solve 
this optimal control problem. Since the control variable, the 
substrate feed-rate occurs linearly in the system, the optimal 
policy is bang-bang and/or follows a singular trajectory. The 
control variable is maintained at its upper or lower bound, in 
bang-bang control, while along the singular arc the control 
variable varies continuously between the lower and upper bound. 

Cazzador (1988) has studied the optimal feed rate profiles 
for fed batch fermentation that maximises the biomass production : 
and the effect of initial conditions. |, 

Fishman and Biryukov (1974) considered the problem of optimal ; 
glucose feed strategy for pencilin production. However, they did 
not give a numerical scheme, and ignored volume changes arising 
from the addition of the feed. 

Weigand et al . (1981) reduced the problem of optimal f eed i 

rate determination for single cell protein production to one in' 
which the total time of operation is minimised for a fixed rate of I 
biomass production. This was analysed for the repeated fed-batch | 
operation. 



Park and Rogers have experimentally studied the factors which 
influence the growth kinetics of the L-phenylalanine strain of 
E-coli NST74. From the batch culture data it was observed that 
the increased initial glucose levels and increased L-phenylalanine 
product concentration caused a decrease in the specific rates of 
growth and rate of product formation, thereby reducing the yields. 

In fed-batch culture a relatively high L-phenylalanine 
concentration was achieved. 

The application of the continuous maximum principle though 
elegant posseses many numerical problems. For problems linear in 
the control variable, the control action is divided into several 
parts - when the variable is at its upper bound, when it is at its 
lower bound or when it varies continuously. The convergence of 
the numerical scheme usually needs a priori information about the 
exact sequence of the different control actions. This is usually 
obtained using information from the physics of the problem (Modak 
(1986) ) . j 

In many problems however the interactions between the ; 
different variables can be very complex. Under these conditions 
it is not possible to decide on the exact sequence of control ' 
profiles which results in optimal performance. It is also , 
experimentally difficult to maintain an optimal profile which ; 
varies continuously. In bio-chemical systems this problem is i 
often avoided by adding the substrate as pulses at discrete 
instants of time. 

Sargent et al . (1984) used a Sequential Quadratic Programming 

approach to determine the optimum feeding policy of a reactant in j 
a series-parallel reaction sequence. They divided the entire j 



duration into different stages and obtained the flow rates 
required for optimum performance in each stage. 

Biegler et al. (1984) have developed a technique based on 
Sequential Quadratic Programming and orthogonal collocation which 
enables us to evaluate the optimal variable efficiently. 

Rosen and Luus (1991) have discussed the methods for 
calculating gradients for nonlinear optimal control problems for 
the rapid convergence. 

In this work we discuss the approximation of the optimal 
policy by discretising the substrate- feed rate with pulses. The 
substrate is assumed to be added at discrete points of time. The 
reaction is assumed to have a Haldane dependency on the substrate 
concentration. When the initial conditions of the reactor are 
stoichiometric, the order of the equations can be reduced by one 
and this permits an analytical solution. For the case when the 
initial state is non-stoichiometric the system of equations has to 
be solved simultaneously. Here we have used a ' Sequential 

Quadratic Programming technique to obtain the pulse feed-rate 
policy for different number of sub- intervals . 

The effect of the initial state of the reactor is determined 
for both these cases. The optimal operation is discussed for (i) 
the final volume is fixed, final time is fixed, (ii) the final 
volume is free, final time is fixed, (iii) the final volume is 
fixed, final time is free, (iv) the final volume is free, final 


time is free. 



CHAPTER 2 


Problem Formulation and Method of Solution 


In our mode of operation the substrate is added 
instantaneously at discrete points of time. These time instants 
signal the beginning of each sub- interval . The entire cycle of 
reactor operation here is divided into different sub- intervals . 
In each sub-interval the reactor is operated in batch mode. The 
pulses which are added at the beginning of a sub-interval change 
the initial conditions for that sub-interval. 

The entire time of operation [t Q ,tg] is divided into N 
different sub-intervals. The evolution in the i sub- interval 
(t^_ 1 < t < t^) is governed by, the equations of the batch reactor 


X 

S 


U(S) .X 



1(a) 

1(b) 


The equations represent the reactor state before the addition of 

the pulse at t^. Here X and S represent the concentration of cell 

mass and substrate respectively and V is the volume of the reactor 
til 

in the 1 interval and i±{S) is specific rate of growth. 

These evolution of the system is subject to initial condition. 
These conditions specify the state at the beginning i.e. t~, or 
before the addition of the pulse at t as 


X'V - X o 

2(a) 

S(t o) 

- S o 

2 (b) 

v(t;> 

“ V o 

2 (c) 


For the numerical work a substrate inhibition model is taken from 
the literature. The dependency of the specific growth rate on S 



is assumed to be given by 


U(S) = 




m 


- Km S 

1 + S- + K. 
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From the mass balance equations at the time instants when the 
pulses are added, we have, 

V, 

I r- 

■i-r 


x(t! <tT_ 1 ) 


3 (a) 


3 + V v Vl s F 


3 (b) 


i i 

where the subscript F stands for feed conditions. Here t^, tT 

represent the instants immediately after and just before the 

addition of the pulse at t = t^. The pulse is assumed to contain 


only substrate and no biomass 


Our objective is to seek the 


pulses A\A's to be added at t^'s which will maximise the biomass 
concentration X(t^). A schematic indicating the different 
variables is depicted in Fig. 1. . ; 


Methods of Solution: 

We solve this problem using two techniques 

(i) Integral method 

(ii) Sequential Quadratic Programming (SQP) 


(i) Integral Method: This technique can be used only when the 

initial conditions are in stoichiometric. The initial conditions 

are said to be stoichiometric when they satisfy 
X(t j 

— + S (t Q ) = S p 4(a) 

After the addition of the first pulse (AV^ the volume becomes 

V n = V + AV. 5 (a) 

1 o l 

The biomass and substrate concentrations after the pulse addition, 



i . e . , at t + , 
o 


are 


X(t o>-^ x «g 

+ V o v - V 

s(t o> - vr 8 <t-> + -V-2 s F 

i. 1 


5(b) 

5(c) 


Adding 1 the two equations 5 (b _ c) we obtain 

r ■*». \ l. i -i 

= 'V^ [ T~ + s (t o } ] 


x(t+) 

— r- + s (t-) 


+* v o r x(t o } 0 . -.1 v o 

Y + S (t o } I + S F " S F 


= S T 


This last identity follows from 4 (a) . 

The evolution of the system (la-b) during the batch is such that 


d f X J\ 

HE [ Y + S j = 0 


or + S(t) = S F 


for t . . < t < t . 
l-l l 


6 (a) 


6 (b) 


clearly 


X(t') ‘ X(tT .) 

— Y — + s (t i ) = s p = — y" + s(t7_ ] 


6 (c) 


Hence the variables X, S satisfy the stoichiometric relation 

for all 't' when the initial conditions at t~ are stoichiometric. 

This can be used to eliminate S in favour of X as 

X 


S = g (X) = S 


F Y 


(7) 


Then equation 1(a) and 2(b) can be reduced to the single equation 


X = n (g (X) ) X 

Integrating this from t* to t 1# we obtain 


( 8 ) 


X(t”) 


I 


dX 


= t„ - t 


Xfi (g (X) ) 1 o 


9 (a) 


X(t 


Thus, the total final time of operation t f we have 



t f fc N " fc o 


X(t") 

n -L 


^ 4- 

X(t^) 


X(t-) 


x (t+) 


dX 

xjTTxl 


X(t-) 


4 * ... + 


X(t N-l> 


9(b) 


Using the jump discontinuity conditions 3(a,b) which arise due to 
the instantaneous pulses this yields 


. V. 
X(t + ) 

t V 1 


V 

X(t 2> Vi 

n' ^ 


x(t o> v: 


X(t*) 




The variables which are to be determined are (i) the volumes 
of the reactor, V^, in each sub- interval i or alternatively the 
pulses to be added at t^'s and (ii) the duration of each 
sub-interval. Alternatively we can evaluate the cell mass 
concentration at the end of each sub-interval, X(t^), instead of 
the duration of each sub- interval 9 s~- ince these uniquely determine; 
the duration of each sub-interval. In the first part of the work 
we consider our manipulated variables to be X(tt) and V^'s 

Our objective is to maximise the biomass concentration at the 
end of the cycle i.e., X(t^) or X^ . The subscript f is used to 
denote the final state of the reactor. The total number of. 
variables to be determined are 2N. These are N, At^s, and N, 


AV i ' s 


At the optimum conditions 


SX(tt) 


for i = 0 


11 (a) 


and 


0 


11(b) 



8 


using these on (10) yields 
4(X(t?)) = Mjx(tt) 


for i = 1 , 2 , ... . (N-l) 


4 


rx v i 

o o 


V, 




Combining both (12a-b) 


i+l> 

for i = 1,2, (N-2) 


4 


fX V 1 
o o 


V, 


- M ( x(t l>) ■ u ( x<t i» - «'( x < t hi» ^f) 


for Isis N-l 


Writing it out for each, i 


4(X(t*)) = u(X(t*)) = ... = d(X(t^_ 1 ) = 4(X(t+)) 

V« n V_ ^ , V, 


4 ( X(t l ) v 1 


= 4 


X ^2^ V ) ~ *'* “ ^ ( X ^N-l^ V 

2 ' v N-l' 


= 4(X(t^)) 


the possible solutions are 


f X V 

x(t+) = x(t+) = ... = xu*^) = -2-0 


x(t*) ^ = x(t^) 


x v 

o o 

V. 


Since the form of kinetic equation is such that there 
utmost two values of X, X^ (* X 2 ) for which 
4 (X 1 ) = 4(X 2 ) , 


A possible solution is 


V, 


V. 


V, 


V, 


V, 


V. , i 

n: i-l 


and X(t+) = X(t 2 ) 


r 

X (t + . ) 

N. “I 


u 


12 (a) 

12 (b) 

12 (c) 

13 (a) 

13 (b) 

14(a) 

14 (b) 

can be 

15(a) 

15 (b) 


This ensures that we have N sub- intervals i.e. no sub- interval is 
of "zero" duration 


t 


f 


V 


1 


ru ru X f 


f di ) . + r 

dX 

, T dx 

J x M(X) + J 

X o V o 

X MX) + 

■ * J X ju(x) 

u 

V 1 



I-’ M 
r • V f 




15 (c) 


15(d) 


ru 

■ J 

X V 
o o 

1-- w 

r V, 


dX 


ru 


X, 


X n(X) 


+ (::-2 


) r 

J Xu 


dX 

(xT + 


r dx 

J ttu 


(X) 


15(e) 


u 


u 


Here the variables to be found are u, r. So we have for optimal 
conditions 

ax- ax f 

= 0 and — z— = 0 


3u 


3r 


This yields 

U(u) = n (ru) 

, X V X 

M (ru) = u \ — 

^ r ^ V f > 

combining (I6a-b) , since we do not 
of zero duration 



The three equations 15(e), 16(c), 

for a fixed V^, t^ to obtain 


16 (a) 
16(b) 

want the first-interval to be 

16 (c) 

16 (d)] 

16 (d) are solved simultaneously 
r, u, Xjc under the optimal 


conditions . 
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(ii) Sequential Quadratic Programming (SQP): 

The integral method proposed can not be applied when the 
initial conditions are non-stoichiometric . This case can be 
analysed using sequential quadratic programming. We discuss the 
solution methodology for two cases 

(i) When the duration is divided N equal sub- intervals . 

(ii) When the duration is divided of each sub-interval is not 
constrained to be equal. 

(i) Equal Intervals: 

Here in fed-batch mode of operation, we divide the duration 

of the operation into n equal time intervals. The substrate feed 

is added at the beginning of each sub-interval. The volume of the 

reactor is constant in each interval. The optimisation problem is 

to determine the optimal feed pulses to maximise final cell mass 

concentration (X^) . We assume cell mass growth to be limited by 

one substrate as done earlier. Here the duration of each 

smte 

sub-interval is known . the total cycle time is fixed. 

A 

The differential equations which govern the system evolution 
are (la-b) . Constraints are imposed on the final reactor volume 
at the final time (t^) , 

V(t f ) = V f (17) 

The updating conditions at the beginning of each interval are the 
conditions (3a-b) . 

In this approach the different charges a to be found for 

the optimal performance are chosen to satisfy the constraint. The 
numerical method used here is Sequential Quadratic Programming 
(SQP) . The convergence of this method needs the accurate 
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determination of derivatives of the objective function with 
respective to manipulated variables, ip this case pulses (AV^) . 

The success of this method depends on the accurate 
computation of these derivatives. We have used state trajectory 
sensitivity method proposed by Rosen and Luus, 1991 to calculate 
the derivatives for faster convergence. 


Calculation of sensitivity with respect to AV^s: 

For the case of equal sub-intervals we recall that 
our objective function is X f and that the manipulated variables are 
AVI ' s for i = 1 . . . .N 

Differentiating la-b with respect to AV_. , we have for j = l,...,N 


ax 

ax 

ax 

+ M 

as 

3AV . 
1 

ax 

3AV. 

1 

as 

3AV. 

1 

as 

as 

ax 

as 
+ as 

as 

3AV . 
1 

ax 

3AV . 

1 

3AV . 

1 

d f 

ax ) 

= 4(S) 

ax 

x ( 4 (s) ) 2 r K m 1 I as 

Ht ( 

3AV. J 

3AV . 

1 

4 m L s 2 K i J , aAv j 


similarly 


d 

r as ] 

1 -n(S) 

ax 

X (4 (S) ) 2 

r -K , 

as 

Ht 

Uav j j 

1 Y 

aAVj 

Y 4 

m 

L s 2 K i . 

3AV . 
1 


j = l,...n,t. _ < t < t. 18(b) 

J 1-1 1 


The four equations (la-b, 18a-b) are to be solved simultaneously 

from t* „ to t. where j = l...n. The following updating 

1-1 f 

conditions 

for j - 1....N are valid for the sensitivity functions 


ax 




3AV. 




Ipi 

v 2 
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as 


v. 


3AV . “ * (S F~ S(t j-i ) ~^~ 2 ~ 

3 J V . 

when j > 1 we should integrate only two equations l(a-b) from 

to t j-i 

for j = 1 . . . .N-l 
i = j+1, N 


3X ,^ + 


V, - V. 


3AV . ' l-l 


(t, j = x(t: , ) 


V. 


i i-i . i-i ax 


l-l' v 2 

i 


V i 3AV. (t i-l } 


as , (s(t i-i> - V (v i' v i-i> 

SAVj = ^2 


V. 

+ ) 
V ± 3AV_. v i-l J 


+ ax(t.) 

Integrating these equations from t to t. T yields the — 

3 N 5AV . 


This will give the gradients of objective function with 
respective to manipulated variables, (the pulses) . 


(ii) Unequal Intervals: 

In this approximation we divide the duration into N different 
intervals. Here we have 2N variables to be found for optimal 
performance. Our manipulated variables are the AV^'s and At^'s. 
For this case we have to find the derivatives with respect to 
AV^'s as well as At^'s. The calculation of gradients with respect 
to AV^'s are the same as discussed above. The calculation of 
gradients with respect to At^'s as follows. 

Differentiating la-b with respect to AV^ , we have 

SX 3X 3X 3X as 

3At j = ax 3At j 3S 3At j 


3S _ 3S 3X 

3At . = 3X 3At . 
3 3 


as 


as 


4 * 


3S 3At j 
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3X 

TTT = JU (S) 

ax 

x X(U( S) ) 2 r K m 1 ' 

as 


0At . ^ V ' 

3 

3 At . 

D 


4 m L s 2 K i - 

aAt . 

J 3 

~ 2 0(a) 

similarly 






d f as 1 

-U(S) 

ax 

X(n(S) ) 2 f” K m 

i 

as 

dt [sAt j J 

Y 

aAt . 

3 

Yju o2 

m L S 

+ IT 
i - 

3At . 

3 




j = 1, . . .n, tj 

< t < 

tj 20(b) 

The four equations 

1 (a- 

■b) , 20 (a-b) , are 

to be 

simultaneously 

integrated from 

t + to 
o 

t f . 

We can do it < 

ctherwise eliminating 


unnecessary computations. 

Suppose if we differentiate with respect to j th time interval 


then integrate the equations from tt to t f . Because 


ax 


dX 




aAt . 

3 

t=t 7 

3 

" dt 

t=t 

j 


as 


dS 




aAt . 

3 


dt 

t=t 

j 


updating conditions 

are 


ax 


V j-1 

ax 

(t j-l 

aAt . 

3 

V. 

3 

3At . 

3 

as 


V j-1 

as 

(f- - 

aAt . 

3 

V. 

3 

aAt . 

3 

(t j-l 


The effect of initial conditions on the optimal feed profile 


and the objective function is also tabulated. 
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Fig . 1 . The variables in the different stages of a cycle. 
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CHAPTER 3 


Result and Discussions 

In this section we discuss the results of the computations 
obtained using the methods discussed earlier. We divide this 
section into two parts. Part 1, consists of results obtained 
using the integral method, non- stoichiometric equal and unequal 
intervals for the case of fixed t f and fixed V f . Part 2 deals 
with the results obtained for remaining three modes of reactor 
operations . 

In our approach we have described the entire period of 
operation into 'N' sub- intervals . The volume of the reactor is 
increased from V q to V f by adding pulses AV t at the beginning of 

t" V» 

the i sub- interval . For the case of equal time intervals the 
problem reduces to determining the N-variables, AV^s to maximise 
X f . 

For the case of unequal sub- intervals, we have to obtain 2N 
variables, the pulses to be added at the beginning of the 
sub- intervals, as well as their durations. 


PART 1: 

The integral method can be used only when the biomass 'X' and 
the substrate concentration 'S' satisfy the stoichiometric 
relationship. This enables us to reduce the two dimensional 
system to a one -dimensional system. The optimal profile of the 
pulses is obtained by the simultaneous.... solution of the three 
equations 15(e), 16(c), 16(d) for u, r and X f . The parameters 

chosen for simulation are K m = 0.397458, = 62.9, 
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a = 0.347695, Y - 0.2, Sp = 10, X = 


P = 10, X q = l.o, V q = l.o, V f = 2.5, 


T f " 5 - 


In Table- 1 we depict the numerical results of the 
computations for different values of 'N' , the number of 
sub-intervals. The first (N-l) sub-intervals are of equal 
duration while the last sub-interval is distinct from these. We 
have also depicted the corresponding u - the biomass 


concentration, r - the ratio of the reactor volumes of 


successive 


intervals i.e., V. _/v., (i 

i+l' i' v 

the end of the duration. 


1, . . . ,n-l) , and the biomass at 


The biomass concentration X^ increases with the number of 
pulses. This confirms our expectations as increasing N results in 
obtaining a better approximation. In the batch mode the biomass 
concentration is 1.7537. 

The improvement in the performance when operated in the 
fed-batch mode in this case is marginal. For the choice of 
parameters listed here the fed-batch reactor gives only a 0.542% 
improvement over the batch reactor for N=10 . Fig. 2a depicts the 
variation of reactor volume for N=5 . The corresponding variation 
of X is shown in Fig. 2b. The addition of the pulses results in 
discontinuities in the biomass concentrations. 


The integral method proposed here has some inherent 
limitations. The method can converge to infeasible solutions.. 
Thus although r > 1, u > 0 and X f > 0, V 1 = V f r 1-n , may be less 
than V , This solution is not feasible since it yields V < V 

O xu 

Such solutions do satisfy the equations but however are physically 


meaningless . 
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A tacit assumption in the integral method is that we add 
pulses at the beginning of each sub- interval . Clearly a solution 
to the equations is that some of the pulses are zero, i.e. , r = 1, 
for some cases, while for others r > l. we have not considered 
this case. This method also implies that the duration of a 
sub- interval is zero when r = l. It does not consider the 
possibility of reactor operation when the i^^ 1 pulse can be zero, 
and the duration of the i th sub-interval is non-zero. 

The optimum value reported as in Table-1 is hence subject to 
these tacit assumptions. 

Sequential Quadratic Programming: 

For the case when the initial conditions do not satisfy the 
stoichiometric relationship we can not use the integral method 
discussed above since now the problem is two dimensional. We have 
solved this problem using Sequential Quadratic Programming. The 

, COCIA 

function 'CONSTR' available in MATLAB used to solve this 

A 

problem . 

Equal Sub- intervals: 

In this method all sub- intervals are considered to be of 
equal durations. We consider the operation where the final 
reactor volume and the duration of the entire cycle are both 
externally fixed. The duration of each sub-interval is hence 
determined once we fix 'N' . Hence in this problem one has to 
determine only the N-pulses, i.e., AV^s which maximises X(t f ) . 

In this mode of reactor operation the optimal profile will be 
dependent on the initial state of the reactor. There are three 
characteristic quantities which characterize the substrate 
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concentration, S Q , the initial substrate concentration, S„ the 

* F 
feed concentration, S the concentration where v(S) has a maximum. 

The relative ordering of the different quantities gives rise 

to different cases. These have been tabulated in Table 2. The 

performance of the reactor i.e., the sequence of the pulses for 

each of these cases is tabulated in Table 3 for fixed V f = 10, 

t £ = 30 and N=5 . 
f 

We now discuss the trends of the optimal profile obtained in 
Table 3, for the different conditions. 

★ 

Case I: S Q < S < S p : The initial concentration is lower than the 

★ . * 

S . Since the reaction rate is maximum at S , the performance is 

optimal when the substrate concentration is maintained closer to 
★ 

S . It is hence necessary to add the first pulse to increase S 

o 

Tk 

to S . The optimal pulse distribution which satisfies all 
constraints is shown in Table 3. The optimal X f here is 2.8372. 
For this case the batch operation has an X f = 1.6012. Fig. 3a 
depicts the volume of the reactor after the addition of each 

pulse. The volume changes dis continuously at the point of addition 
of each pulse. The corresponding biomass and substrate 

concentrations are shown in Fig. 3b(i) and Fig. 3b(ii). 

•Jc "k 

Case II: S p < S < S Q : In case 2a, S Q is approximately S and S p 

is very low. Addition of S p will hence hurt the reactor 

performance. The pulse addition is hence delayed till the end. 

This addition is necessary since we have to satisfy the volume 

★ 

constraint. For case 2b, S Q is significantly higher than S . 

Addition of S,, reduces to S* now and this yields the optimal 
profile shown , Here the fed-batch performance is same as batch. 
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Case III: S Q < S F < S : Since the initial concentration is low 

addition of the feed increases the substrate concentration close 
to S*. Here it is advantageous to add all the pulses in the first 
sub- intervals as this brings S Q closest to S* . In this case the 
final biomass concentration is same for fed-batch and batch 
reactors . 

Case IV: S < S p < S q : Here the first pulse is added as in case 
III, at the first sub-interval. This brings S q closest to S* and 
the reaction rate is closest to the maximum. Here again the 
performance of the fed-batch and batch reactors are identical. 

★ 

Case V; S < < Sp : Here the initial concentration is closer to 

★ 

S than Sp. Addition of the feed would thus increase S and lower 
c o 

the reaction rate. This addition should be hence delayed as far 
as possible. This confirms the trend shown in the table. 

★ 

Case VI: S p < S Q < S : Here again addition of the feed would 

drive the concentration away from S*, lowering the reaction rate. 
Hence this should be delayed. In this case the substrate is added 
in the last interval to satisfy the volume constraint . 

Unequal Sub- intervals (Fixed t f and Fixed V f ): 

In Table-4 we depict the results of the computations for the 
case of unequal time intervals. Here t^ = 30, and = 10. We 
have presented the results for the same N as in Table-3. This 
enable us to compare the two results. 

The qualitative trends of the two Tables agree well. The 
performance of each case for a given 'N' is better when the 
Interval u u to of unequal duration as compared to when w<: 
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restrained the intervals to be equal. 

We simulated the behaviour of the system for the case of 
stoichiometric initial conditions with SQP algorithm for unequal 
sub-intervals. The parameters set chosen is the one used to 
generate the results using the integral method. The SQP method 
also yields the same optimal profile as the integral method. The 
length of (N-l) sub-intervals is equal and the last one is 
distinct. The successive volumes V i+1 /V. also bear a constant 
ratio for i-l,2,...,N. The optimal objective function also 

matches to within a desired numerical accuracy with the results in 
(Table-1) . 
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PART 2 : 

We now discuss the application of this method to remaining 
three modes of reactor operation (i) fixed final time, free final 

free final time and fixed final volume (iii) free 
final time, and free final volume. 

The optimisation problem solved so far have some cases where 
the volume of the reactor is made up in the last sub- interval . 
This is forced on the system since it has to satisfy the 
constraint on the reaction volume. It would hence appear that 
should this constraint be removed the optimal profile be one in 
which no pulses are added. To test this hypothesis removed the 
constraint on the reaction volume. We now discuss the results for 
the situation when the reactor is operated over a time t f and the 
final free volume . 

For cases 2, 4, 6 of Table-5, here, there is no pulse added 
for optimal performance. The optimal performance is in batch mode 
with reactor volume equal to the initial volume. hence the final 
pulse added for the fixed volume problems, in Table-4 is only to 
satisfy the constraint of the fixed final volume. 

For cases 1, 5 pulses are added in each sub- intervals . The 
final volume for both the cases is less than 10. In case 1: = 
7.96, while for case 5 it is 7.1925. The corresponding cases in 
Table-4, also have non- zero pulses being added in each 
sub- interval . 

For case 3, the optimal performance is found when the initial 
pulse is non-zero as in Table-3 For this case the optimisation is 
carried out by constraining the final volume V f to be lower than 
the upper bound. The upper bound chosen for this case study wusi 

central library 
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20. The reactor is filled in the first pulse, since this raises 
★ 

S Q closer to S . This is again identical to the optimal profile 
in Table-4 . 


One would expect a similar trend for case 4, i.e., addition 
in the first pulse to being S Q to S . However, for the values 
chosen S Q , S p are very closer to S* and hence batch mode at V q is 
opt 1 mal . 

The improvement of the performance as measured by the 
concentrations X f for each case in Table-5 is much better than for 
the corresponding cases in Table-4. The free volume problem has a 
better performance than the reactor with a fixed final volume. 

The optimal performance for a given number of pulses may not 
be necessarily obtained when the sub- intervals are of the same 
duration. The constraint of each sub-intervals results in 
obtaining a sub-optimal profile for a given N. 

To obtain a time optimal profile for a fixed number of 
intervals N, we now discuss the results obtained when the 
sub- intervals are of unequal durations. Now for a given number of 
N sub-intervals we have to determine 2N variables, the N At^'s and 
the N AV^'s. We solve this problem also using "CONSTR' of MATLAB. 
The convergence of the algorithm is based on the accurate 
calculation of the sensitivities of the objective function with 
respect to the control variables. These are computed as discussed 
earlier . 

In Table-6, we depict the results for fixed t f and free V f , 
with unequal sub-intervals. The effect of initial conditions on 
the optimal profile is similar to that predicted in Table-5. For 
each case the optimal solution is better when we relax the 
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const 3r8.int s of equa.1 sub- intervals as expected. 

Table-7, contains the results for V f = 10, and free t f . The 
solution method requires an upper bound on free time. This was 
fixed at 150. In most of the cases as shown the optimal profile 
was found such that the cycle time was 150. A characteristic 
feature of most of the runs in that the system posses multiple 
optima . 

The convergence to the different optima depends on the 

initial guesses supplied. Although the profile can be very 
different there is no change in the objective function. For this 
mode the improvement in the objective function over the batch mode 
is negligible in most cases. This occurs because the system tends 
to choose the largest possible time to allow the reaction to go to 
completion. This results in maximum biomass. 

Table -8 depicts the results for the free V f and free t f 

problem. Here again we have to impose upper bounds on these 

variables to the optimal profile obtained in most cases is when 

the reactor is operated in the batch mode. The choice of t £ , V f 

here is such that the biomass X £ at the end of the cycle is much 

more than that obtained in the earlier modes of operation for each 

* 

combination of S q , S , S p . The corresponding biomass and 

substrate concentrations for unequal intervals are shown in 
Fig.4b(i) and Fig.4b(ii). 

Table-9 contains the results for V £ = 10, and free t f for 
equal intervals of time. We impose an upperbound on free time. 
This was fixed at 2500. The convergence to the different optima 
depends on the initial guesses supplied. Here also the 
performance of fed-batch reactor is almost equal to the batch 
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mode . 

Table-10 depicts the results for free V f and free t 
Here we impose an upper bound on these variables. The 
profile obtained in most cases is when the 
the batch mode. 


problem. 

optimal 


reactor is operated in 
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CHAPTER 4 

Conclusions 

Here we have demonstrated different methods which can be used 
to obtain the optimal performance of fed-batch reactions. The 
flow rate over the entire cycle is approximated by discretised 
pulses . 

The integral method is applicable only when the initial 
conditions are stoichiometric. It is also tacitly assumes that 
all pulses are non-zero. 

Sequential quadratic programming allows us to obtain optimal 
profiles for the case of non-stoichiometric initial conditions. 
The effect of S , Sp, S has been studied and explained 
qualitatively. 

The numerical convergence of SQP algorithm requires the 
accurate computation of the sensitivities. An efficient way to 
evaluate these is discussed. 

For a particular set of parameters the case of stoichiometric 
unequal intervals yields the best possible optimal solution. The 
SQP method for the stoichiometric case can yield optimal profile 
even when the integral method gives infeasible solution. 
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Table 1 


Integral Method 


N 

(N-l) At' s 

At 

n 

u 

r 

X f 

batch 

(N=l) 

0 

5.0 

- 

- 

1.7537 

3 

1 . 1818 

2.6363 

0 . 8125 

1.4252 

1.7623 

5 

0 . 6714 

2.3142 

0.8952 

1.2231 

1.763 

7 

0.4675 

2.1954 

0.9277 

1.1505 

1.7632 

10 

0.3207 

2.1136 

0.9508 

1.1010 

1.7632 
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Table £ 


Parameters used for Table ,3-10) are ^ - 0.208, K m = 2.516, 
K * = 25 ’ Y = °' 8 - X o = »-l. V 0 = 7, S* e 7 . 93095 


Case 


0 . 142857x10' 


2a 

8 

0.14 

2b 

20 

4 

3 

0 . 142857X10" 3 

5.0 

4a 

8 

7.95 

4b 

20 

10 

5a 

7.95 

75 

5b 

10 

75 

6 

3.0 

1.0 


75 


Description 


S o < S ■= S F 

S < s* < S 

F O 

F O 

S o < S F < S * 
< S F < S o 
S * < S F < S o 

S * < S o < S p 

S* < s n < s 
o F 

s p < S o < S* 
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Table 3 


Non-stoichion&etr ic , Equal Intervals of Time, Fixed T^., Fixed V^. 


Case 

AV 1 

AV 2 

AV 3 

av 4 

av 5 

Fed-batch 

X f 

Batch 

X f 

1 

0 . 8343 

0 . 0226 

0 . 0487 

0 . 1057 

1 . 9888 

2 . 8372 

1.6012 

2a 

- 



- 

3.0 

2.7845 

2.5542 

2b 

3.0 

- 

- 

- 

- 

2.4624 

2.4624 

3 

3.0 

- 


- 

- 

0 . 5493 

0.5493 

4a 

1.5294 

1.4706 

- 

- 


3.0911 

3 . 0911 

4b 

3 . 0 

- 

- 



2.212 

2.212 

5a 

0.0047 

0.00224 

0 . 0486 

0.1064 

2.818 

2 . 6703 

1.1918 

5b 

- 

- 

- 

- 

3.0 

2.5857 

1.11945 


6 


3 . 0 


1.157 


1.001 
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Table 4 



Non- 

stoichiometric. 

Unequal Intervals 

Fixed T f 

and Fixed V f 





i 





Case 

1 

2 

3 

4 

Fed-batch 

5 x f 

Batch 

X f 

1 

AV. 

l 

0.8577 

0 . 1122 

0.1518 

0 . 1847 

1.6926 

3 . 1817 

1.6012 


At i 

15 . 805 

6 .692 

4.31 

3.20 

- 



2a 

AV. 

- 

- 

- 

- 

3 

2 . 967 

2 . 5542 


At. 

l 

8 . 1178 

8.1195 

8.1195 

5.6432 

- 



2b 

AV. 

X 

3.0 

- 

- 

- 

- 

2.4624 

2.4624 


At^ 

5.578 

5.9645 

6.095 

6.158 

6.204 



3 

AV i 

3.0 

- 

- 

- 

- 

0.5493 

0.5493 


At i 

5 . 849 

5.9744 

5.975 

6.036 

6.170 



4a 

AV. 

0.762 

0.761 

0 . 76 

0 . 717 

- 

3 . 0911 

3.0911 


At . 

l 

1.968 

3 . 968 

• 5.968 

7.979 

10.117 



4b 

AV. 

l 

3.0 

f - 

- 

- 

- 

2.2120 

2.2120 


At. 

l 

7.64 

7.78 

7.78 

6.81 

- 



5a 

AV x 

0.028 

0 . 1131 

0.1518 

0.1851 

2.522 




At i 

15.7815 

6.697 

4.315 

3.21 

- 

3 . 1808 

1.198 

5b 

AV. 

l 

- 

- 

- 

0.2271 

2.7729 

3.1248 

1.1195 


At i 

9.354 

8.804 

7.447 

4.3956 




6 

4Vi 

- 

- 

- 

- 

3.0 




4t i 

9.881 

9.9806 

10.081 

0,0587 

-• 

1.17895 

1.0096 
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Table 5 


Non-stoichiometric, Equal Intervals of Time, Fixed T^., Free V^. 


Case 

AV 1 

AV 2 

AV 3 

AV 

4 

AV 5 

Fed-batch 

X f 

Batch 

X f 

1 

0.8343 

0.0226 

0.0487 

0 . 0953 

0 

3 . 92424 

3.8938 

2X 

- 

- 

- 

- 

- 

4.23861 

4.23861 

2 b 

- 

- 

- 

- 

- 

2.6805 

2.6805 

3 

10.094 

- 

- 

- 

- 

0.8469 

0.8469 

40- 

- 

- 

- 

- 

- 

4.23861 

4.23861 

4b 

- 

- 

- 

- 


2 . 6805 

2 . 6805 

5a 

0.0046 

0.0226 

0.0487 

0.1060 

0.0106 

4.36098 

4.3077 

5b 

- 

- 

- 

- 

_ 

4.4159 

4.4159 


6 


4.23861 4.23861 



Table 6 


Non-stoichiometric, Unequal Intervals Fixed T f and Free V f 


Case 


1 

2 

3 

4 

Fed 

5 

-batch Batch 

X f X f 

1 

AV i 

0.8419 

0.0432 

0.0547 

0.0644 

- 

3.9245 

3 . 8924 


At. 

l 

10.312 

5.4711 

3.7455 

5.3416 

5.1308 



2a 

AV. 

X 

- 

- 

- 

- 





At . 

l 

5 . 806 

5 . 9048 

6.0029 

6.0976 

6.1876 

4.2386 

4.2386 

2b 

AV. 

l 

- 

- 



- 




At i 

5.9989 

5.9990 

6.001 

6.001 

6.002 

2.6805 

2.6805 

3 

AV i 

10.094 


- 


- 

0.84696 

0.84696 


At i 

2.5047 

6.0501 

6.8642 

7.2006 

7.3804 



4a 

AV. 

X 

- 

- 

- 

- 

- 

4.2386 

4.2386 


At i 

5.801 - 

5 . 901 

6.001 

6.101 

6.1964 



4b 

AV. 

- 

- 

- 

- 

- 




At. 

1 

5.4717 

5 . 8039 

6.1004 

6.2575 

6.366 

2 . 6805 

2.6805 

5a 

AV i 

0.0134 

0 . 0482 

0.0614 

0.0727 

- 

4.36719 

4.3 057 


At . 

l 

10.8834 

5.6831 

3.8212 

3.7499 

5.8624 


j 

5b 

AV i 

- 

- 

- 

- 

- 




At i 

7.673 

7.673 

6.226 

4.9555 

4.067 

4.4159 

4.4159 

6 

AV ± 

- 

- 

- 

- 

- 




At. 

6.003 

6 . 003 

6.003 

6.003 

5.988 

1.68421 

1.68421 


Table - 7 


Unequal Intervals of Time, Fixed Final Volume, and Free Final Time 


Case 




i 



Fed-batch 

Xj 

Batch 

Xf 


1 

2 

3 

4 

5 

1 

AVi 

0.2601 

2.7399 




18.07 

18.07 


Ati 

10.119 

9.503 

21.064 

9.7196 

10.689 



2a 

A Vi 


«. 


0.8429 

2.1571 

4.55 

4.55 


AU 

30.016 

30.116 

30.216 

30.189 

29.462 



2b 

AVi 

1.4175 

0.6373 


- 

0.9452 

12.23 

12.23 


A t{ 

28.827 

29.09 

29.34 

30.26 

32.49 



3 

AVi 

3.0 

- 


_ 

- 

1.2708 

1.2708 


AU 

- 

2.0399 

3.91 

52.72 

91.34 



4a 

AVi 

1.5854 

0.8037 

0.204 

0.204 

0.204 

6.458 

6.458 


AU 

29.8 

29.9 

30.0 

30.1 

30.2 



4b 

AVi 

1.1024 

0.2938 

* 

0.2059 

1.3979 

13.67 

13.67 


AU 

29.511 

29.356 

29.554 

30.2 

,31.73 



5a 

AVi 

0.7125 

2.2875 


- 

- , 

22.522 

22.522 


AU 

29.8195 

23.05 

3.953 


1.264 



5b 

AVi 

3.0 

- 

_ 

_ 

- 

23.67 

23.67 


AU 

80.88 

41.36 

1.27 

" 

11.83 



6 

AVi 

. 1.6178 

1.3822 


_ 

- 

1.99 

1.99 


AU 

45.813 

38.49 

31.15 

23.29 

11.254 





Table - 8 


Unequal Intervals of Time, Free Final Volume, and Free Final Time 


i Fed-batch Batch 

Case 1 2 3 4 5 Xf X f 


1 

A Fi 

0.0154 

1.8217 

0.0113 


Af{ 

2.5 

48.361 

5.281 

2a 

AVi 

- 

_ 

- 


A U 

16.37 

16.47 

16.57 

2b 

AVi 


_ 



A h 

14.12 

14.22 

14.32 

3 

AVi 

146.685 

0.6511 

4.6065 


AU 

161.36 

162.77 

163.77 

4a 

AVi 

- 

- 

- 


AU 

13.461 

13.56 

13.66 

4b 

AK 


_ 

_ 


AU 

•15.45 

16.45 

17.45 

5a 

AVi 

0.6945 

3.884 

0.1977 


AU 

29.0819 

42.09 

20.74 

5b 

AV 

0.694 

3.88 

0.197 


AU 

29.1 

42.1 

20.7 

6 

AVi 

_ 

_ 

- 


AU 

45.0 

38.0 

31.0 


1.2379 

6.953 

89.8138 

70.744 

55.807 

11.76 

16.67 

16.77 

6.5 

6.5 

14.42 

14.52 

16.1 

16.1 

14.098 

164.77 

26.8596 

166.1 

3.8635 

3.8635 

13.77 

13.86 

8.1 

8.1 

18.45 

19.43 

16.1 

16.1 

0.7982 

10.693 

87.4255 

41.407 

56.2522 

22.36 

0.798 

10.69 

. 87.43 
41.0 

56.37 

21.96 

23.0 

11.0 

2.5 

2.5 



Table - 9 


Fqual Intervals of Time, Fixed Finn! Volume, and Free Final Time 


Case 

A t 

AVi 

av 2 

A^3 

W A 

av 5 

Fed-batch 

Batch 

Xf 

1 

56.45 

3.0 

- 

- 

- 

- 

18.0701 

18.0701 

2a 

45.954 

- 

- 


0.52 

2.48 

4.55 

4.55 

2b 

17.098 

1.8931 

0.9501 

0.1567 

- 

- 

12.23 

12.23 

3 

169.927 

3.0 

- 

- 

- 

- 

1.2701 

1.2701 

4a 

192.236 

0.0001 

0.0004 

0.0003 

0.0003 

2.9988 

6.458 

6.458 

4b 

14.684 

- 

0.4368 

- 

- 

2.5632 

13.67 

13.67 

5a 

67.811 

- 

1.7762 

1.2238 

- 

- 

22.52 

22.52 

5b 

55.474 

3.0 

- 

- 

- 

- 

23.67 

23.67 

6 

44.691 

0.186 

1.8868 

0.927 

_ 

_ 

1.99 

1.99 


Table - 10 


Equal Intervals of Time, Free Final Volume, and Free Final Time 


Case 

A t 

A Vi 

av 2 

av 3 

av 4 

av 5 

Fed-batch 

X/ 

Batch 

x f 

1 

30.0 

4.82 

- 

5.41 

9.4 

73.4 

55.81 

52.27 

2a 

16.76 


- 

- 

- 

- 

6.5 

6.5 

2b 

29,1 


- 

- 


- 

16.1 

16.1 

3 

29.97 

93.0 

- 

- 

- 

- 

3.727 

3.7065 

4a 

16.06 

- 

- 

- 


- 

6.5 

6.5 

4b 

15.89 

- 

- 

- 

- 

- 

16.1 

16.1 

5a 

29.77 

42.022 

0.004 

0.005 

1.536 

49.43 

56.252 

33.416 

5b 

29.948 

0.0663 

2.1801 

- 

8.9635 

81.7902 

56.367 

35.561 

6 

24.24 



- 

_ 

- 

2.5 

2.5 



•• 
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Fig. 3a : Volume profile for the fixed V^ and t^ with equal 


intervals . 
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Biomass concentration profile for the fixed V f and 
tj with equal intervals. 


Fig. 3b(i) 



Substrate concentration profile for the fixed V f 


Fig. 3b(ii) : 


with equal intervals. 
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Fig. 4b(i) 


t^ with unequal 





Fig. 4b(ii) : 


Substrate concentration profile for the free 
and tj with unequal intervals. 



